Quantum correlations of two optical fields close to electromagnetically induced 

transparency 
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We show that three-level atoms excited by two cavity modes in a A configuration close to elec- 
tromagnetically induced transparency can produce strongly squeezed bright beams or correlated 
beams which can be used for quantum non demolition measurements. The input intensity is the 
experimental "knob" for tuning the system into a squeezer or a quantum non demolition device. 
The quantum correlations become ideal at a critical point characterized by the appearance of a 
switching behavior in the mean fields intensities. Our predictions, based on a realistic fully quan- 
tum 3-level model including cavity losses and spontaneous emission, allow direct comparison with 
future experiments. 

PACS numbers: 42.50.Dv,42.50.Gy,42.65.Pc 



Using optical nonlinearities it is possible to manipulate 
optical beams to the level of quantum fluctuations, pro- 
ducing squeezed states [l| that are important resources 
for quantum information with continuous variables 
Related to the field of squeezing is that of quantum 
non demolition (QND) measurements on an optical field, 
vifhere quantum correlations between two different modes 
of the electromagnetic field are exploited to overcome the 
back-action noise of a quantum measurement (sf . Besides 
the fundamental interest in the theory of measurement, it 
was shown that QND correlations of propagating beams 
have direct application in quantum communication pro- 
tocols as teleportation The best single back-action- 
evading measurement on optical beams was performed 
using cold atoms inside a doubly resonant cavity Q . We 
suggest that these performances could be significantly 
improved by tuning the system close to the electromag- 
netically induced transparency (EIT) conditions 

Already in the nineties, theoretical studies showed that 
a lambda three-level medium close to EIT conditions in 
a cavity can be used to obtain squeezing [7]. Contrarily 
to previous proposals, here we assume that two different 
modes are resonant in the cavity. For small and sym- 
metrical detunings from the upper level of A three-level 
atoms (see FiglT]), absorption is suppressed and the dis- 
persive non linear response gives rise to a rich scenario 
where either self correlations (squeezing) or cross QND 
correlations can be established in the output beams. The 
correlations become ideal at a critical point that we char- 
acterize analytically. The technique we propose is exper- 
imentally accessible, and first experimental steps in this 
directions were done in atomic vapors without a cavity 
. Here we show that the presence of the cavity is a cru- 
cial advantage especially if one can reach the good cavity 
limit. 

We consider N atoms in a cavity excited by two modes 
symmetrically detuned from the upper level of a A three- 
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FIG. 1: Two cavity modes interact with the atoms in a A 
configuration close to EIT conditions. 



atomic transition. We define A, — ^ — ^ the atomic de- 
tunings normalized to the decay rate of the optical coher- 
ences 7^, = (71 -|- 72)72 where 71 -I- 72 is the total popula- 
tion decay rate of the upper level; 9j = '^"^^ '^^ the cavity 
detunings normalized to the cavity decay rates Kj, and 

N 

Cj — ^ the cooperativities where gj are the coupling 
constants for the two considered transitions. We use nor- 
malized variables proportional to the intracavity and in- 



put fields Xj 



and 



respec- 



tively, where Tj is the (input-output) mirror transmissiv- 
ity for the field j . We name v and w the normalized polar- 
izations between levels 1-3 and 2-3 v = —[\/2/N){R^), 
w = —{^/2/N){S~) where R and S are collective oper- 
ators constructed from the single atom operators |1)(3| 
and 1 2) (3 1 as in 9]. The master equation and the semi- 
classical equations describing the A system with two cav- 
ity fields, with the same notations introduced here, are 
given and discussed in detail in [lo| where this model was 
successful to reproduce the experimental results of 

Let us consider a set of parameters symmetric for the 



two transitions: 



\yj\ 



\v\ 



7i 



7, Kj 



9j = (empty cavity resonance for both fields), and let 
Ai = ~A2 = e be small and positive. In Fig. [2] we show 
in rescaled units the stationary intensities of the intra- 



level scheme as in FigdJ For j = 1, 2 let ojj be the fre- cavity fields Ij 



■/4Ce as a function of the common 



quency of field j andhtOaj the energy of the corresponding intensity of the input fields Y = \y\ /ACe. With a solid 



2 




FIG. 2; Stationary intensities of the intracavity fields Ii, and 
I2 as a function of the common intensity of the input fields 
Yi = Y2 — Y . In solid line the solution Ji = h- The thick 
(red thin) line correspond to stable (unstable) solutions. In 
dashed-dotted line one of the two stable solutions with 7i 7^ 
l2- Parameters: 2e = 0.125 and Ci = C2 = 250, 71 = 72 = 
IOki, K2 — k.1, 9i — 62 = 0. 



line we have plotted an S'-shaped solution with Ii ^ l2- 
A stable branch of this solution appears for y > 1. The 
negative slope branch and the lower branch very close to 
zero intensity are both instable and play no role in the 
following. For Y < 1, apart from the solution Ii = I2, 
we get two other solutions with Ji l2- In the figure 
we show one of them with Ii > l2- The second one is 
obtained by exchanging Ii and 12- Both solutions are 
stable in the considered case 6*1 = 6*2 = 0. 

We choose now two values of the input intensity, in 
turn above and below the turning point Y = 1, and 
show the stationary solutions for intracavity fields inten- 
sities as the cavity detunings vary in Fig. [3l The sta- 
ble branches of these curves (thick lines) can be easily 
obtained experimentally by sweeping the cavity length 
[lo| . We vary 9i and 62 keeping them always equal which 
would imply the use of two driving fields of close optical 
frequencies AA/A <ti 1 (and for example different polar- 
izations). For Y = 1.05 i.e. 5% above the turning point 
(upper half of Fig. [3]) the stable solutions for the intra- 
cavity intensities are Lorentzian-looking curves symmet- 
rically shifted by a small amount from their empty-cavity 
positions for both fields. Only for 61 — 62 — the two 
fields have the same stationary amplitude in the cavity 
corresponding to the stable high-transmission branch of 
the S'-shaped curve in Fig[2l For Y — 0.95 i.e. 5% below 
the turning point (lower half of Fig. [3]) the situation is 
rather different: the stable solution for the two fields 
switches between a high-intensity and a low-intensity 
curve being always Ii ^ I2 although |yi| — \y2\- In con- 
trast with the previous case this situation is very far from 
the independent-fields FIT solution and the fields arc in 
fact strongly coupled. 

Let us now introduce the usefuU correlations to carac- 
terize the quantum fluctuation properties of the system. 





FIG. 3: Intracavity field intensities Ii (left half), I2 (right 
half) across the cavity scan. Upper half: Y = 1.05. Lower 
half: Y = 0.95. The thick (red thin) lines correspond to stable 
(unstable) solutions. The other parameters are as in Fig. [2l 



il^e**^, the squeezing spectrum is defined as 



Sfiu) = 1 + 2k, / e"'^'{: SxUt) 5xUo) : ) dt (1) 



where SXl^ denotes the time dependent fluctuation of 
the operator X"^ around a steady state point. The col- 
umn indicates normal and time ordering for the product 
inside the mean. Sj^ = 1 is the shot noise and 5*^^ = 
means total suppression of fluctuations in the quadrature 
Xj*". The crossed correlations between the two fields are 
described by the coefficients Cs, Cm and Vs\rn till] char- 
acterizing a QND measure of the amplitude quadrature 
-^in q£ Qj-^g field, the signal, performing a direct measure- 
ment on the phase quadrature y°"* of the other field, 
the meter. Among the three coefficients Cg quantifies 
the non-destructive character of the measurement, 
its accuracy and refers to the to the "quantum state 
preparation" capabilities of the system. 

Cs = C(X'",X°"*) , = C(X™,y°"*), (2) 

Vs\m = (i-c(x°"*,r°"*)) (3) 

where for two operators A and B we define 

\{A,B)\^ 



C{A,B) = 



/+00 
-00 



{A,A){B,B) 



with 



(4) 



'-{A{t)B + BA[t))dt. (5) 



For a given quadrature of the j field: X 



ajC 



Superscripts in and out refer to the input and output 
cavity fields. By calling (/i™ and 0°"* the phases of the 
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input and output fields in steady state, and choosing 
field 1 as the meter and field 2 as the signal, we define 

QND measurement Cm = Cs = 1, and Vs\m = 0. 

The quantum fluctuations counterpart of Fig. [3] (top) is 
shown in Fig. [4] (top) where squeezing of the output fields 
optimized with respect to the quadrature S'j'^** (w = 0) is 
plotted as a function of the cavity detuning. A large 
amount of squeezing is present in both fields close to 
di = 0. As one can see from Fig. |3] (top) the two fields are 
well transmitted by the cavity for 6i = 0, and the system 
efficiently converts the input coherent beams into bright 
squeezed beams. Correspondingly to Fig. [3] (bottom) for 
Y — 0.95, in Fig|4] (bottom) we plot the coefficients Cg, 
Cm and Vs\m across the cavity scan. The useful quan- 
tum correlations are calculated by a linearized treatment 
of quantum fluctuations around the stable stationary so- 
lution as in [10] . Despite the fact that the two fields have 
different intracavity intensities at 9i — 0, they play here 
symmetrical roles for the QND scheme; the figure corre- 
sponding to the reversed scheme 1 2 being obtained 
by reflection of the plots 6i ^ —Oi. 




FIG. 4: Top: Best squeezing of the fields across the cavity 
scan for Y = 1.05 and uj = 0. Squeezing of field 1 (2) is 
plotted with a solid (dashed) line. Bottom: QND coefficients 
across the cavity scan for Y = 0.95 and u> = 0. Cm (red 
dashed-dotted line), Ca (solid line), K|m (blue dashed line). 
Parameters as in Fig[3] 

We show in Fig. [5] the frequency dependence of the 
quantum correlations both below and above the turning 
point y = 1, for a fixed value of the cavity detuning close 
to zero. For values of the cooperativity parameters cur- 
rently obtained in experiments, QND coefficients such as 




FIG. 5: Top: Squeezing spectra of field 1, for Y — 1.05 in 
the center of the cavity scan {9i = 0.0013). Best squeezing in 

solid line and amplitude squeezing S-^ ^ in purple dashed line. 
Bottom: QND spectra for Y = 0.95 and Oi = 0.0018. C™ (red 
dashed-dotted line), Cs (solid line), Vsi^ (blue dashed line). 
The other parameters are as in Fig[2] 



in this regime, representing a significant improvement 
with respect to previously obtained resul ts p based on 
the so called "ghost transition" scheme [l3 |T [loj . Al- 
though we concentrate here on the good cavity limit, in 
which as we will show the quantum correlations become 
ideal approaching the turining point Y = 1, some QND 
correlations between the two modes persist also in the 
bad cavity limit. For example for e = 0.25, C = 25, 
K = 37, 6*1 = 6 X 10-^, Y = 0.9 and uj = O.I7 we get 



Cs — C„ 



0.72, Vs 



0.26. 



In the limit of weak atomic detunings, useful analytical 
results can be obtained. The analytical solution of the 
semiclassical equations of the system at steady state is 
given in [l3| • By expanding the steady state polarizations 
V and w to the first order in e we obtain 



V = i 



4:exi\x2\ 



(|a;i|2 + |x2|2)2 



w = —i- 



Aex2\xi\ 



\x2\^r 



(6) 



By inserting ^ in the equations for the intracavity fields 
amplitudes, with = \y\, 9j = and Cj = C, we 
obtain at steady state a "universal solution" for rescaled 
field intensities. For Y < 1 there are two stable solutions 



Y 



(i±^); h 



Y 



(iTry) 



(7) 



where rj — \/l — Y^. For Y > 1, out of two solutions 



/ ; / 



1± Wl 



Cs = 0.98, C„ 



0.95, K 



0.05 can be acheived 



Y2 



(8) 
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the one with the plus sign is stable and the other one 
unstable. Solutions ([Zl)-® are indistinguishable from 
those of the full three-level model in Fig. O The 
phases of the input fields with respect to intracavity 
fields (which are taken real at steady state) are 0™ = 

for y < 1, and 



at an 



— — atan 



= atan[l/2/] = -0^" for Y > I. For the output 



fields 



lout 
"2 



/j™ in both cases. 



In order to study the quantum properties of the system 
analytically we further assume that (i) 7 3> k so that the 
atomic fluctuations follow adiabatically the field fluctu- 
ations, and (ii) the noise from spontaneous emission is 
negligible, which we found true when the cooperativity 
is large enough. In this limit, using the steady state po- 
larizations ([6]), we can solve analytically the equations 
for the fleld fluctuations and obtain the correlation func- 
tions. 

For Y > 1 and Ii = I2 = I and taking k^^ as the unit 
of time, we obtain 



Sxi 



-6x-i 



21 ^ 



J = 1,2. 



(9) 



These equations describe two independent two-photon 
processes for which instabilities and squeezing have been 
studied extensively [l^]- The best squeezing spectrum 
for each field is 



4a 



S'^^'ico) = 1 2 



yielding perfect squeezing at zero frequency at the turn- 
ing point where y = 1, / = 0.5 and a = 1. 

For Y <1 and h 7^ I2 the fiuctuations of the two fields 
are coupled. For Ji > I2 we get 



5X1 
5Y^ 



-SXi 



1 



-6Y1 + i- 



1 



Y 
+ rj 



5Y^ 



Y 



5Xi - 2i TJSX2 



(11) 
(12) 



The equations for field 2 are obtained from (|lip and ([72]) 
by changing the sign in front of 77 and of i. Simple ana- 
lytical expressions can be obtained for the squeezing and 
the conditional variance V^\.„ of the fields at a; = 



Qtnt 



gbest ^ ^ 



Qphase 

^3 



-3+4 (13) 
(14) 



showing that the fields have diverging phase noise and 
become perfectly correlated at the turning point. We 
checked that the spectra in Fig. [5] are well reproduced 
by the analytical results. 

In conclusion, in a symmetrically detuned FIT scheme, 
and for equal input intensities Y of the two fields we have 
shown the existence of a universal ^-shaped steady state 



curve (Fig[2]) which divides the parameter space into two 
parts: for input intensities higher than the upper turn- 
ing point of the curve, the quantum fluctuations of the 
fields become quadrature dependent and can be reduced 
in a quadrature, while for input intensities lower than 
the turning point, crossed phase- intensity quantum cor- 
relations build up between the two fields. The system 
becomes a perfect "squeezer" or an ideal QND device at 
the turning point. The "universal" point Y = I, can 
be identified experimentally by the appearance of the 
switching behavior described in Figl^l and can be used 
as a reference in the parameter space to choose either the 
squeezing or the QND effect and to optimize it. An im- 
plementation using either a vapor ISjl, or a trapped cold 
atoms in an optical cavity [5|.|l4|.|15| seems within the 
reach of present technology. 
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